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ABSTRACT 

We discuss observations of the magnetic field, column density, and turbulence in 
the Cold Neutral Medium (CNM). The observed quantities are only indirectly related 
to the intrinsic astronomical ones. We relate the observed and intrinsic quantities by 
relating their univariate and bivariate probability distribution functions (pdfs). We find 
that observations of the line-of-sight component of magnetic field do not constrain the 
pdf of the total field B to t very well, but do constrain the median value of B to t- 

In the CNM, we find a well-defined median magnetic field 6.0 ± 1.8 [iG. The CNM 
magnetic field dominates thermal motions. Turbulence and magnetism are in approxi- 
mate equipartition. We find the probability distribution of column density N±(HI) 
in the sheets closely follows N±(HI)~ 1 over a range of two orders of magnitude, 
0.026 < N±(HI) < 2.6 x 10 20 cm~ 2 . The bivariate distributions are not well enough 
determined to constrain structural models of CNM sheets. 
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1. INTRODUCTION 

Beginning February 1999 we used the Arecibo telescope 1 to begin a series of Zeeman-splitting 
measurements of the 21-cm line in absorption against continuum radio sources. Heiles & Troland 
(2004; Paper III and references therein) discussed technical details of the observing technique 
and data reduction and presented observational results, including magnetic field measurements. 
As described in Paper III, our analysis of the data identified Gaussian components in the HI 
absorption spectra. We assume that each component samples a single, isothermal, sheet-like region 
of CNM. For 69 such components, we are able to estimate (1) the line-of-sight magnetic field 
strength Bi os , often subject to significant error, (2) the line-of-sight column density N(HI)i os , 
and (3) the contribution of turbulence Vturb los to the line-of-sight velocity dispersion. Systematic 
instrumental errors are small and the uncertainties are Gaussian distributed. Therefore, our survey 
has yielded a statistically well-defined ensemble of observed values for each of these three CNM 
quantities. 

From these data, we seek to answer several astrophysically significant questions. For example, 



lr The Arecibo Observatory is part of the National Astronomy and Ionosphere Center, which is operated by Cornell 
University under a cooperative agreement with the National Science Foundation. 
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what are the probability distribution functions (pelfs) of magnetic field strengths, column densities 
and turbulent energies in the CNM? Are any of these quantities statistically related to each other? 
If so, can we determine if the magnetic field lies preferentially in the planes of CNM sheets or 
perpendicular to them? What is the ratio of magnetic to turbulent energy in the CNM? As we will 
see below, the data provide answers to some, but not all, of these questions. 

We observe quantities that are only indirectly related to the intrinsic astronomical ones. For 
the magnetic field, we observe the line of sight component Bi os , not the total field B tot . For the 
column density N(HI), we must account for the fact that the Cold Neutral Medium is in sheets, 
not spheres; our observed column densities Ni os are always larger than the intrinsic column density 
normal to the sheet N± . To properly interpret our results and to address the questions listed above, 
we must consider how the intrinsic and observed quantities are related in a statistical sample. These 
statistical transformations, and the results of applying them to our data, are the focus of the current 
paper. 

We treat both univariate distributions and bivariate distributions. The univariate distributions 
are most interesting because we obtain definitive results; in contrast, the possible correlations that 
could be revealed by bivariate distributions are obscured by noise and inadequate numbers of data. 
Accordingly, the reader who is interested in astrophysical results can concentrate on the sections 
dealing with univariate distributions of magnetic field, column density, and turbulent velocity. We 
recommend beginning with §2, which introduces the notation and concepts of intrinsic and observed 
quantities. We develop the theoretical relationships between observed and intrinsic univariate 
distributions in §3, and use these to obtain the actual intrinsic univariate distributions in §3.1. 
Finally, we discuss astrophysical implications in §7, §8, §9, and §10. 

2. SOME IMPORTANT WORDS ON NOTATION 

We need to introduce important distinguishing nomenclature because of the anisotropic nature 
of magnetic-field-related quantities. First, we make the crucial distinction between line-of-sight 
(los) quantities, such as Bi os , and the intrinsic astronomical ones such as Btot- We will often refer 
to the latter with one the shorter terms "intrinsic" or "astronomical". The important point is to 
distinguish between observed and intrinsic quantities. 

The line-of-sight component is the quantity to which the telescope responds, and we designate 
these with the subscript los. There is an additional complication for the magnetic field, which is 
the presence of significant instrumental noise. Thus, we must distinguish between the line-of-sight 
component of the field, designated by Bi os , and the actual measured value, designated by B b s . We 
define 



Bios = the line— of —sight component of B tot ; 



(la) 
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B bs = the observed field strength . 



(lb) 



The essential difference between Bi os and B Q b s is 



Bobs — Bios + SB) 



noise 



(lc) 



where 5B no i se is the uncertainty contributed by random measurement noise. 

Noise is small enough to neglect for column density and velocity, so the observed quantities are 
essentially identical to the los ones; we will use the subscript los for the observed quantities, with 
the implicit assumption that their uncertainties from noise are negligible. For column density, we 
always refer to HI so we often write N instead of N(HI). We usually consider sheets, for which the 
column density perpendicular to the face is N(HI)±_ and for which the apparent column density 
along the line of sight is 



where 9 is the angle between the sheet's normal vector and the line of sight. Similarly, for velocity 
linewidths, we use the symbol Vt ur b,ios to indicate the line-of-sight component of the turbulent 
velocity. 

Finally, the probability density functions (pdf) of intrinsic quantities, such as B tot , differ from 
those of the line-of-sight quantities. We will always use the symbol 4> f° r the pdfs of intrinsic 
quantities, and the symbol ip for line-of-sight (los) or the observed ones. Unless otherwise specified, 
units are always as follows: magnetic field B in //Gauss; column density N in HI atoms 10 20 cm -2 ; 
velocity V in km s _1 . 

3. THEORETICAL CONVERSIONS OF INTRINSIC ASTRONOMICAL 
PROBABILITY DENSITY FUNCTIONS TO OBSERVED ONES 

We begin with some light theory by considering elementary transformations of magnetic field 
and column density in statistical distributions. First, we consider how the intrinsic probability 
density function (pdf) of Btoti defined as 4>(Bt t), converts to the observed histogram or pdf of 
observed B[ os , defined as ip(Bio S ), under the assumption that fields are randomly oriented with 
respect to the observer. Note the important distinction: Btot is the total field strength; Bi os is the 
observed line-of-sight component, which is always smaller. 

Next we incorporate one of the fundamental results of Heiles & Troland (2003; Paper II), 
namely that the CNM components are thin sheets. We define N± = N(HI)± as the HI column 
density perpendicular to the sheet; the observed quantity is Ni os = N(HI)i os , which is always larger. 
As with the magnetic field, we consider how the pdf (j>(N±) converts to the observed histogram 
or pdf ip(Ni os ), again under the assumption that the sheets are randomly oriented with respect 



N los = N(HI) los = N(HIUs 



N(HI) ± 



(2) 



cos 9 
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to the observer. We also consider the statistical transformation of the distribution of the actual 
nonthermal velocity dispersion Vt ur b to that of the observed one I4ur&,Zos under the assumption that 
turbulence is perpendicular to the magnetic field. 

We then tackle the bivariate distributions. First we consider how the bivariate distribution 
<P(Btot, N±) converts to the observed one i/j(Bi 0S , Ni os ). We then assume two extreme models, one 
with the fields always perpendicular to the sheets and one with fields parallel to the sheets, and 
assuming random orientations derive the observed bivariate distributions ip(Bi os , Ni os ) for the two 
cases. We illustrate and discuss the result by considering delta-function distributions of B tot and 
N±, and we also apply the transformation of the observed tp(Bi os , Ni os ) to its intrinsic counterpart 
<p(Btot, N±). Finally, we examine the bivariate distributions involving the pairs {V tur b,ios-> Ni os ) and 
(Vturb, los, Bios), which produces little in the way of useful results. 



3.1. Conversion of the univariate distributions 

3.1.1. Conversion of the intrinsic d>(B to t) to the observed ip(Bi os ) 

We first consider the simple case in which all clouds have the same B tot , which is randomly 
oriented with respect to the observer. The line of sight component Bi os is 



Bios = Btot cos 9 , (3) 

where 6 is the angle between the field direction and the line of sight. 6 can run from to ir, but 
the intervals from to ^ and | to ir are identical except for a change of sign in Bi os . It's simpler 
and no less general to consider the smaller interval 6 from to | so that we can ignore the slight 
complications of the sign change. In this case, the pdf of 6 is the familiar 



<j>e(9) = sm6 



(4) 



and we wish to know the pdf of Bi os , which is given by (see Trumpler and Weaver 1953 for a 
discussion of these conversions) 



Vb{Bi os ) = MQ(Bios)} 



d[6{B los )} 
dBios 



(5) 



which gives 




I n if < Bios < Btot , c s 

ip(Bios) = { * tot . (6) 

otherwise 
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In other words, B[ os is uniformly distributed between the maximum possible extremes and B to t 
(actually ±Bt t)- This leads to the well-known results that in a large statistical sample for which a 
constant Btot is viewed at random angles, both the median and the mean observed field strengths 
are half the total field strength and that Bf os = — jp. More generally, observed fields are always 
smaller than the actual total fields, and with significant probability they range all the way down 
to zero. 

Now suppose B^t has an arbitrary pdf (p(B to t)- The univariate pdf (f)(6) becomes the bivariate 
pdf (f)(Btot,&), and we assume B to t is independent of the observer's location so that (p(B to t,0) = 
4 > B tot (Btot)4>8(0)- Note that we introduce subscripts on the different 0's to distinguish them, instead 
of designating them with different Greek letters. To obtain ip(Bi os ) we again follow the standard 
techniques; it's easy to integrate over 9 and obtain 

mios)=r ^IdBtot, (?) 

where the symbol [Bi os > Btotmin] means the larger of the two quantities. The presence of Btot i n 
the denominator means that smaller ranges of B[ os are emphasized. This is commensurate with the 
equation 6's uniform pdf for a single field value. We note that this can be regarded as an integral 
equation for (p(B to t) and it is straightforward to invert. 

Figure 1 illustrates the solution of equation 7 for four functional forms of <p(B to t) plotted against 
[jfy^p where the subscript 1/2 denotes the median value. These forms include the following: 

1. <p a Kronecker delta function (DELTA FCN), (pB tot {Btot) = S(B to t - B to tfl), yielding tp a flat 
function (as discussed immediately above, equation 6); 

2. <f> a flat distribution (FLAT FCN) between < \B tot \ < B , yielding tp oc In (j^J ; 

3. (p a weighted Gaussian (EXP FCN), 



yielding ip a Gaussian with dispersion Bq. 

4. <fi a Gaussian (GAUSS FCN) with dispersion Bq, yielding tp oc E\ (^-^tj, where E\ is the 
exponential integral of order 1. 

All four (p(Btot) are plotted with respect to B Bt °l /2 ' so ^ e m edians of all lie at unity on the 
x-axis. However, the means differ. Similarly, the medians and means of the associated i/j[Bi os ^ 
differ from each other. These relationships between median and mean are summarized in Table 1. 
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DELTA FCN 
FLAT FCN 
EXP FCN 
GAUSS FCN 




E T0t'/ ID T0T,1/2' 



DELTA FCN 
FLAT FCN 
EXP ECN 
GAUSS FCN 



D LOS'' ID T0T,1/S' 



Fig. 1. — Top panel: The intrinsic cj)(B tot ) for four representative functional forms. Bottom panel: 
their line-of-sight counterparts ip{Bi os ). The vertical scales are arbitrary. 

The medians and means for ip{Bi os ) are all about half those for (Bt t), which is a direct result of 
the weighting by B^ ] in equation 7. 

Figure 1 is disappointing from the observer's standpoint, because the observed distributions 
ip{Bi os ) do not differ very much. These differences become smaller — inconsequential, in fact — when 
one includes measurement noise, as we shall discuss in detail for our data in §4.1. Unfortunately, 
given the inevitable errors in any observation that is sensitive to B[ os , it seems practically impossible 
to distinguish among different functional forms for (f>(Btot)- Nevertheless, the average value of Bi os 
is close to half the average value of B tot for a wide range of intrinsic pdfs of the latter; this also 
applies to the medians, but less accurately. Therefore, this rule of thumb may be used to estimate 
the median or average B to t from an ensemble of measurements of Bi os . 



3.1.2. Conversion of the intrinsic <j>(N±) to the observed ip(Ni os ) for sheets 



Here we assume the CNM is distributed in sheets having HI column density N± in the direction 
perpendicular to the sheet. If the normal vector to the sheet is oriented at angle 9 with respect to 
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the line of sight, we have 

Nios = (9) 
cos 

To find tp(Ni os ) we follow the same procedures as above in §3.1.1. For a single value of N± we 
obtain 



f if N los > N± 

mios) = ^ r (io) 

[ otherwise 

For a single N±, iVj os has a long tail extending to infinity. The median value of Ni os is N loSjl / 2 = 
2Nj_, reflecting the increased observed column for tilted sheets. The mean value of Ni os ((Ni os )) is 
not defined because the integral diverges logarithmically; of course, this doesn't occur in the real 
world, where sheets don't extend to infinity. For example, if all sheets have aspect ratio 5:1, then 
(N los ) = 1.6N ± and N losA/2 = 1.7N ± . 

For an arbitrary pdf (j)(N±) and infinite slabs, we obtain 

maxl 



*P(N los ) = ^2" / N ± <f>(N ± ) dN ± . (11) 



As above, this can be regarded as an integral equation for (f)(N±) and it is almost as straightforward 
to invert. 



3.1.3. Conversion of the intrinsic <j)(Vturb) t° the observed tp(Vt ur b,ios) for turbulence 

perpendicular to B to t 

The observed velocity width comes from two sources, thermal and nonthermal. We can separate 
these because we have independent measurements of the kinetic temperature. Thus, we can derive 

Table 1. Medians and Means of B tot and Bi os for Representative pdfs 





Btot, 1/2 


(Btot) 


Bl os ,l/2 


(Bios) 


DELTA FCN 


1.00 


1.00 


0.50 


0.50 


FLAT FCN 


1.00 


1.00 


0.40 


0.52 


GAUSS FCN 


1.00 


1.18 


0.38 


0.59 


EXP FCN 


1.00 


1.04 


0.44 


0.51 
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the line-of-sight (los) nonthermal ("turbulent") line width and its associated energy density. How 
this relates to the total turbulent energy density depends on whether the turbulence is Id, 2d, or 
3d. 

Here we assume that this turbulence is 2d, i.e. we assume that it is restricted to motions per- 
pendicular to the mean magnetic field. We make this assumption because we show below that the 
typical turbulent Mach number is equal to 3.7. For velocities that are parallel to the magnetic field, 
such turbulence would produce strong shocks that would damp very rapidly. However, velocities 
that are perpendicular to the magnetic field can be as high as the Alfven velocity without pro- 
ducing shocks; this is the basis for considering turbulent motions perpendicular to the mean field. 
We caution, however, that numerical simulations of magnetohydro dynamical turbulence find that 
magnetic fields do not ameliorate turbulent dissipation (MacLow et al. 1998), so our assumption 
might not have any basis in physical reality; in this case the turbulence is isotropic (3d) and the 
observed distribution i/j(Vt urb j os ) is equal to the intrinsic one (f)(Vt urb ). 

Suppose that the magnetic field is oriented at angle 9 with respect to the line of sight, as 
in §3.1.1. Suppose that turbulent motions are perpendicular to the fieldlines and isotropic in the 
azimuthal directions around the fieldlines, and along one direction perpendicular to B have width 
Vturb', this makes the full turbulent width 2 1 / 2 Vt UT b- Then line-of-sight width is 



Vturb,ios — V turb sm9 . (12) 
To find tp (Vturb, los) we again follow the same procedures as in §3.1.1 and obtain 



turb, los) 



V 2 - V 2 

turb turb,los 



-1/2 



if Vturb,los < Vturb 

otherwise 



(13) 



Given a value for V turb , ip(V turbilos ) 
well-defined with 



oo as Vt U rb V tur b los j but the cumulative distribution is 



cum{V turblos ) = 1 



j ( Vturb,los \ ' 
\ Vturb J 



1/2 



(14) 



which gives the median V tur b,ios,i/2 = Q-&7V turb median and mean (Vt ur b,ios) = 0.79V turb . These 
high values reflect the large fraction of sheets tilted to the line of sight, where V tur b,ios is large. 

For an arbitrary pdf 4>{Vturb) we obtain 



^ {Vturb, los) = j ^Tp^ly turb ~ V turb, los] ^ <t>V turb {V turb )dV turb (15) 
JV turb<lo3 v turb 
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In contrast to the two cases above, this integral equation is not straightforward to invert. It is 
a Volterra equation of the first kind and, following the identical example in Trumpler & Weaver 
(1953), it can be rewritten as Abel integral equation. The analytic solution is 

<t>(V t urb) = -^f^-^— T7^~~~ [VLb,los ~ VLb] ^ mturb,los)dVturb,los (16) 

VT dVturb J Vturb Kurb,los 

3.2. Conversion of the intrinsic bivariate distribution (p{B to ti N±) to the observed 

We have measured both Bi os and Ni os and wish to know the bivariate distribution tp{Bi os , N[ os ). 
We proceed by first assuming that B to t, N±, and of course 9 are all uncorrelated. To proceed we 
consider two different models. 

3.2.1. Case of B to t perpendicular to the sheet 

We refer to this as the perpendicular model. The only angle involved is 6, so that the original 
pdf is the trivariate (f)(B tot , Nj_, 0). This case is simplified because only cosO is involved, which 
makes Bi os oc j^—. Converting the original distribution (j)(B to t, N±,9) to the one involving the 
observed parameters yields 

i>(Bios,N los ,N ± ) = ^Btot (^^ ££ ) <f>N ± (M±) (17) 

Here we have chosen to eliminate Btot and express the result in terms of iVj_ ; we could have gone the 
other way. To obtain ip in terms of only the observed quantities, we need to integrate over N±, but 
we cannot do this without knowing (j>N_ L (N±). Later we use the one obtained from observations. 

For now we consider the illustrative case for which all Btot and iVj_ are identical, i.e. (f>B to t {Btot) = 
5 {Btot — Btot,o) and 4>n ± {N±) = 5{N± — N± ), where 6 is the Kronecker delta function. This is a 
trivial case because all observed points fall on the line 

Bios = Btot,oirr^ (18) 

^los 

which is shown in the top panel of Figure 2. 
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Fig. 2. — The theoretical observed joint pdfs ijj(Bi os , iVj os ) for the illustrative case of 5-function 
distributions for Bi os and Ni os . The top panel shows the pdf for Btot perpendicular to the sheets; 
it degenerates into a single line. The bottom panel is for Btot parallel to the sheets; contours are 
spaced by factors of 2 with arbitrary scaling, and the dashed line shows the median Bi os versus 
AW 
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3. 2. 2. Case of B to t parallel to the sheet 

This case is more complicated because Bi os depends on two angles. These are 6, which is the 
tilt of the sheet with respect to the observer, and $, which is the azimuthal angle of the field within 
the sheet. We have Bi os = B tot sin sin <£. This means the pdf <p is quadrivariate. Again we assume 
everything is uncorrelated and eliminate 6 and Obtaining the observed distribution is somewhat 
cumbersome but yields the surprisingly straightforward result 

^(B los ,N los ,B tot ,N ± ) = J^[(B tot N los ) 2 - (B tot N ± ) 2 - (B^N^) 2 }- 1 / 2 ^^)^^) 

(19) 

Later we use observed distributions, but for now we again consider the illustrative case of delta 
functions for B tot and iV|_. Integrating over B to t and N± yields 

^Btos, N los ) = J^[( Btotfi N los ) 2 - (B totfi N ±0 ) 2 - (B los Ni os ) 2 }' l l 2 (20) 
which is illustrated in the bottom panel of Figure 2. 



3. 2. 3. Discussion of Figure 2 

Figure 2 exhibits the joint pdfs for the two sheet models (B tot perpendicular and parallel to 
the sheets — the "perpendicular" and "parallel" models). The median observed column density 
Nlosi/2 i s twice the assumed iV_i_ and the median observed magnetic -B; osl / 2 is half the assumed 
Bt t] these univariate medians are indicated by squares on the two panels. The significance of 
these squares is that half the observed Bi os , and half the observed Ni os , are smaller and half larger. 
Finally, the dashed line in the bottom panel exhibits the median Bi osl / 2 versus Ni os ; we calculate 
this by extracting the conditional pdf i^{Bi os \Ni os ) versus Ni os and calculating the median from its 
cumulative distribution, thus obtaining Bi os l / 2 versus Ni os . 

Figure 2 illustrates a crucial observational signature at large Ni os that distinguishes between 
the two sheet models. More specifically, for the perpendicular model, large Ni os goes with small 
Bi os , and vice-versa for the parallel model. For the perpendicular model, all of the datapoints 
having Ni os above its univariate median (indicated by the square) have B[ os below its univariate 
median. In contrast, for the parallel model most (66%) of the datapoints have Bi os above its 
univariate median: as A^; os gets large, the conditional pdf ip(Bi os \ iVj os ) — > j^b, (B to t 2 — Bi os 2 )~ 1 / 2 , 

los 

which produces the median Bi os i/ 2 — > 0.71B to t- 
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3.3. Conversion of the intrinsic to observed bivariate distributions involving 

turbulent velocity 

In this section, we derive and display the observed bivariate distributions involving ip(Bi os , Vt U rb,ios) 
and tp(Ni os , V tur b,ios), under the assumption that the turbulent velocity is perpendicular to the mag- 
netic field. 



3.3.1. Conversion of the intrinsic 4>(BtotiVturb) to the observed tp(Bi os ,Vt ur b,ios) 



The joint pdf ip(Bi os , Vt ur b,ios) is the same for both the parallel and perpendicular models, and 
indeed does not depend on whether there are sheets or not, because Vturb is perpendicular to Btot so 

21 V2 



that there is a unique relationship = 



( Vturb,los \ ' 
Vturb ) 



. As before, we begin with trivariate 



distributions involving 6, Btot and Vturb- We eliminate 6 to obtain 



1p(Bl os , Vturb. 



r 

,los) = / 



Vturb. 



los 



turb,los 



V 2 

turb 



V 2 

turb,los 

~V 2 

turb 



-1/2 


( 


(PB tot 


Bios 



V 2 

turb,los 

~V^ 

turb 



-1/2N 



0Vt urb {Vturb) dVturb 

(21) 



The top panel of Figure 3 displays this joint pdf for delta-function distributions of Btot and 
Vturb- Because we assume that the B to t and V tur b are perpendicular, strong fields go with small 
velocities and vice-versa. There are no contours; they all collapse into a line because, for the 
delta- function distributions, there is a one-to-one relationship between Bi os and V tur b,ios- 



3.3.2. Conversion of the intrinsic <fi(N±, Vturb) to the observed ip(Ni OS i Vturb,ios) for the 

perpendicular model 

In contrast to the above case, the relation between (p(N±, Vturb) an d ip(Ni osl V tur b ,ios) does 
depend on the model. For the model with B to t perpendicular to the sheet we have Ni os = and 
V t urb,ios = Vturb sin 9. We obtain 



mios,V t 



turb.los) 



N lo 



N± 

N los 



Nj 



-1/2 



V 



turb.los 



Ni 



-1/2N 



cf> N± (N ± )dN ± (22) 



The middle panel of Figure 3 displays this joint pdf for delta-function distributions of -/Vj_ 
and Vt ur b. There are no contours; they all collapse into a line because, for the delta-function 
distributions, there is a one-to-one relationship between Bi os and Vt U rb,los- Because we assume that 
the Btot and V tur b are perpendicular, high HI columns go with small velocities and vice-versa. 
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3.3.3. Conversion of the intrinsic (p(N±, Vt ur b) to the observed cf>(Ni os , Vt ur b,ios) for the 

parallel model 

The parallel model is more complicated, as it was in §3.2.2, because $> enters explicitly: because 
V tur b is assumed to be perpendicular to B to t, V tur b,ios = — (sin # sin (&) 2 ] 1 / 2 . The final 

expression equivalent to 22 is quite complicated, so we deal with it numerically. The bottom panel 
of Figure 3 displays this joint pdf for delta- function distributions of N± and Vturb', the contours 
come from a Monte-Carlo calculation. 



3.3.4- Summary 

The complicated nature of the the univariate distribution ip(Vt ur b,los) m equation 15 means 
that the bivariate distributions that involve Vt ur b,ios are even more complicated and preclude closed- 
form solutions when we use the observationally derived pdfs. However, we can easily present the 
bivariate results for delta-function distributions of B tot , V tur b, an d N±, shown in Figure 3. Two, 
^{Bios-, Vturb, los) an d ip(Ni os , Vturb, los) i degenerate into lines because the observed variables depend 
only on 6. 

4. DERIVATION OF INTRINSIC ASTRONOMICAL UNIVARIATE PDFs 

FROM OBSERVED HISTOGRAMS 

4.1. Derivation of the intrinsic (f>(Btot) from the histogram of observed B a b s 

Figure 4, top panel, exhibits the histogram of measured field strengths B b s - It contains 69 
measurements, of which only 12 have the measurement error b~Bi os ,m < 2.f>Bi os , rn . This histogram 
is not symmetric: it has 42 instances of B Q b s > and 27 instances of B b s < 0. For a randomly 
distributed angle between line-of-sight and the field direction, one expects the numbers to be equal; 
the probability that we obtain this imbalance, or worse, is given by integrating the binomial (coin- 
tossing) pdf and is equal to 0.042. This is fairly low, but is hardly low enough to rule out the 
random distribution. However, the small probability might indicate that the selection of sources in 
the Arecibo sky does, in fact, involve non-randomness. 

We proceed on the assumption that the distribution of angles is, indeed, random. This allows 
us to compare the histogram of the absolute value of measured field |-B; OSim | with the theoretical 
expectation for various assumed intrinsic pdfs of Btot- Figure 4, bottom panel, exhibits this distri- 
bution, but symmetrized so that every entry for a positive B Q b s is matched by one with negative 
B bs- The solid line is best fit the Gaussian pdf 
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Fig. 3. — Bivariate distributions for the model having turbulent velocities perpendicular to the 
magnetic field, discussed in §3.3. These distributions assume ^-function distributions of the astro- 
nomical parameters Vt ur b, Btot, and iVj_. Top panel, the locus of contours of the bivariate distribu- 
tion ip(Bi os , Vturbjos) (observed quantities). Middle and bottom panels: contours of ip(Ni os , Vt ur b,ios) 
Squares show the univariate medians. 
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Fig. 4. — Top panel: the histogram of observed field strengths B Q b s (in fJ,G). Bottom panel: the 
same data, but symmetrized with each observation represented a second time with the opposite 
sign. The smooth curve is the Gaussian fit of equations 23. 
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1 



p obs f^^obs,0^ 



(23a) 



2irBl 



obs,0 



with 



B obsfi = 5.2 ±1.3 a*G . 



(23b) 



To find this best fit and error for B b s ,o, we numerically sampled a range of trial B obs fi values, 
and for each trial value compared the calculated cumulative distribution with that of the data 
by performing the Kolmogorov-Smirnov test. This provides sets of the K-S statistic D and its 
associated probability PKs(B bs,o) that the assumed distribution matches the observed one (see 
Press et al 1997). We determined the best fit for -B o f>s,0 by choosing the one that maximizes PkSi 
which is Pks = 1-0 (meaning that a Gaussian is an excellent fit), and we defined its uncertainty as 
being where the Pks fells to 32% of its peak value (thus mimicking the definition of the la error 
for Gaussian statistics). 

Of course, the 5.2 ± 1.3 fiG of equation 23b represents the pdf of B t, s . This is not the same as 
the pdf of Bios because of observational noise, which is very significant. Therefore, the dispersion 
of Bi os is considerably smaller. Of course, it is Bi os , not B Q bs, which is the quantity of interest, 
so we need to statistically account for the observational noise. This is not straightforward because 
the 69 measurement errors are all different. If they had all been identical, then we could have used 
Gaussian statistics and convolutions to derive the true dispersion Bi os o from the measured B o bs,0- 
Because the errors are not identical, we instead use a Monte Carlo analysis and employ the actual 
measurement uncertainties instead of their rms. 



Here we use the observed values of B Q b s , together with the uncertainties dB no i se , in a Monte 
Carlo simulation. We begin with four different possibilities for the functional form of (f)(B tot ); each 
of these is characterized by the median total field strength Btot, and we calculate 46 uniformly 
spaced field possibilities ranging from 1 to 10 /j,G. 

We assume random orientation with respect to the line of sight. We assume the observational 
errors 5B no i se to be Gaussian-distributed. For each test possibility, we perform the set of 69 
observations many times (20000 "trial runs"). Each time, for each of the 69 measurements we 
randomly generate a value of Btot according to the assumed 4>(Bt t)', orient it randomly with 
respect to the line of sight; derive its observed value Bi os ] generate an uncertainty 5Bi os from a 
Gaussian pdf whose dispersion is equal to the associated value of SB no i se ; and record the resulting 
value of the observed B Q b s = Bi os + SBi os , which includes measurement errors. For each trial run, 



4.1.1. Monte Carlo Method 
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all four possibilities for the functional form use the same values of B tot , orientation, and 5Bi os . We 
make histograms of the absolute values B b s and use the KS test to compare with the cumulative 
histogram of our data. 




Fig. 5. — The KS fits of the Monte Carlo simulations for the four trial functional forms to the data, 
with the median \B tot l / 2 \ (in /j,G) as the independent variable. 

Figure 5 displays the KS fits to the four functional forms. As anticipated in §3.1.1, all four KS 
fits are reasonably good, showing that the data cannot distinguish among these functional forms. 
However, medians for all four functional forms are quite similar, so the data do select the median 
with a reasonably small uncertainty. 

To collapse the results of four functional fits into a single number for the median total field, 
we simply average the fields for the four peaks to obtain 6.0 fiG. To obtain the uncertainty, we use 
the maximum spread of the 32% points, which are 4.2 and 7.8 fiG. These combine to make our 
derived total median field strength 

|5 tot;1/2 | = 6.0±1.8/xG. (24) 

4-1-2. A Single Uncertainty SB no i se for Purposes of Convolution 

Just above, in §4.1.1, we derived the median value B tot i/ 2 = 6.0 /uG. Under the assumption 
that the distribution of observed values without noise is Gaussian, then the pdf of B to t is the 
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weighted exponential of equation 8. In this equation, the median B totl / 2 = 6.0 fiG corresponds to 
B = 3.9 fiG. 

This value of B = 3.9 /iG is really Bi os q, i.e. it refers to what we would see without noise; in 
contrast, with noise we found -Bo.obs = 5.2 ^G. Thus, for this particular pdf for Btot, our ensemble 
of B i, s is equivalent to convolving the ensemble of actual values Bi os with a single Gaussian of 
dispersion 

SB noiseiSngl = (5.2 2 - 3.9 2 ) 1 / 2 = 3.4 ^G . (25) 

Below, in our interpretation of bivariate histograms of our results, we will convolve the Monte- 
Carlo-derived values of Bi os (which have no noise) with a Gaussian of this dispersion to obtain 
simulations of the measured histograms involving B i, s (which includes noise). This is a quick, 
approximate way to illustrate how measurement errors affect the histograms of observed values. 



4.2. Derivation of the intrinsic <fi(N±) from the histogram of observed Ni os 

We first distinguish between two distinct sets of data for Ni os and Vi os . One is the entire sample 
of CNM components from Paper II. The second is the sample we have been discussing, namely that 
comprising the 69 statistically interesting measurements of B b s ; these are drawn from Paper III, 
restricted to those having small uncertainties 5B no i se . This smaller sample has a bias eliminating 
small column densities, which are irrevocably associated with large uncertainties in SB no i se . 

We should not use the smaller dataset of 69 samples for discussing the statistics of column 
density and velocity, because it is biased. Instead, we use the larger dataset, but with some 
restrictions. We include all components from Paper II that satisfy the restrictions that |6| > 10°, 
V? ii > 0, and Vt ur b,ios < 5.25 km s -1 . The first restriction helps to eliminate blending of CNM 
components; the second eliminates 5 components for which the errors happen to yield T^ max < T^; 
and the third eliminates 3 outliers to the analytic approximation below. All this leaves us with 138 
for the sample from Paper II. In the ensuing discussion, we will specify which dataset we are using. 



4-2.1. A Preliminary: latitude dependence of N± 

In §4.2.2 below we will assume that the sheets are oriented randomly and derive the intrinsic 
4>(N±) from the histogram of observed Ni os . Before doing this, however, we consider whether the 
distribution is, in fact, random. In particular, we might expect the sheets to lie parallel to the 
Galactic plane, in which case we would expect Ni os oc ^j^, so we discuss this possibility first. 

Figure 6 exhibits the latitude dependence of Ni os . Looking at the points, one sees some higher 
values near smaller \b\, leading one to suspect that this might be a statistical trend. To test this 
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Fig. 6. — Diamonds are datapoints, iVj os (in 10 20 cm~ 2 ) versus |6|. Each star represents the median 
of Ni os within the 10° wide bin centered on the star; the stars are connected by the dashed line. 
The dotted line is Ni os = si ^j b | , and is meant only to guide the eye. Two points, located at 
(\b\,N(HI) = (28,12) and (19,13), are excluded from the figure (but not the medians) to save 
space. 

suspicion we examine the latitude dependence of binned medians. The errorbars on the starred 
medians were calculated using the absolute values of the residuals, i.e. the error is ~^tjtt^Y)~ > 
where Ni os y 2 is the median and M is the number of points in the bin. The dashed line connects 
the medians of Ni os for 10°-wide bins in \b\. 

If Ni os oc ^jn, then these medians should follow the same dependence; we provide the dotted 
line, which is Ni os = ^jn, for a visual comparison. Visually, there is absolutely no tendency for 
the medians to follow ^jn- The presence of datapoints with high iVj os at small |6| might be real 
or it might be a statistical artifact. If it is real, it might mean that sheets are in fact aligned with 
the Galactic plane to some degree and/or have their own intrinsic distribution of N±. 

There is more to this relationship than random statistics, namely our nonuniform and nonran- 
dom sampling of the sky. Of course, we are restricted to Arecibo's sky, which is a 38°-wide swath in 
declination. In the vicinity of right ascension 04^30™ it passes through the Taurus/Perseus region 
at I ~ 180°, covering a large spread in latitude and a small range in longitude. Large column 
densities persist here up to \b\ ~ 50°, which is an anomaly in the overall latitude distribution of the 
neutral ISM. Thus in Taurus/Perseus, if we were to plot the total observed column density towards 
any direction (as opposed to Ni os , which is the observed column density in one CNM Gaussian 
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component) versus |6|, it would not behave anything like ^jn- On the contrary, in other longitude 
ranges the situation differs. For example, in the vicinity of £ = 240° the ISM column density is 
anomalously small. Thus, plotting quantities versus \b\ also plots them versus a biased distribu- 
tion in I. Accordingly, given the nonuniform sampling of the sky, it is very difficult to establish a 
latitude dependence. 

From the above, we conclude that we cannot extract a latitude dependence of Ni os from the 
data. Accordingly, we proceed under the simplest assumption, namely that the sheets are randomly 
oriented with respect to the observer's line of sight. 



4-2.2. Derivation of 4>(N±) assuming the sheets are randomly oriented 

The top panel of Figure 7 exhibits the histograms of observed column densities Ni os together 
with a curve that, over most of the region, goes roughly as iV^ 1 (see equation 27 below). The 
solid histogram is for the full Paper II sample of 138 datapoints. The eyeball tells us that this is a 
reasonably good fit. Using equation 11, we find that if i^(Ni os ) oc Nr^, then <^(iVj_) has the same 
dependence, namely <f>(N±) oc N^ 1 . Unfortunately, its integral diverges so it is not a valid pdf. 
Rather, we must impose lower and upper limits on this pdf. 

To proceed, we assume 



[ otherwise 

where k = [ln(iV ± max I 'N \min)]~ l an d derive the corresponding ip(Ni os ) from equation 11. This 
yields 





K W 

Los 



Nios < N ± 

min 

N± min < Ni os < N± max 



(27) 



to 



Nios > N ± 



We then numerically cover a grid of trial values for (N± min , N± max ) and for each combination 
derive the corresponding ip(Ni 0S ). From that, we calculate the cumulative distribution and compare 
with the observed one by performing the K-S test, from which we obtain the bivariate probability 
PKs(N± m i n , N± max ) that the assumed distribution matches the observed one. We determine the 
best fit parameters by choosing the combination (N± min , N± max ) that maximizes Prs- 

Figure 8 shows contour plots of ^s^imitii-^imm)' The top panel is for the 138 datapoints 
from Paper II, which is the larger and unbiased sample. The highest contour encloses a well-defined 
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Fig. 7. — Analytic fit and observed distributions for Ni os (in 10 20 cm -2 ). The solid histogram in 
the top panel is for the full set of 138 sheets from Paper II; the dashed is for the restricted set of 
69 from Paper III, used in the present paper. The smooth curve is ip{Ni os ), the pdf of equation 27, 
which is the analytic fit to the 138 datapoints; this curve closely follows Nj~* The bottom panel 
shows the three equivalent cumulative distributions; the diamonds are the solid histogram and the 
squares the dotted one. 
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Fig. 8. — Contours of PKs(N± min , N± max ) (units are 10 20 cm -2 ). The stars are the adopted 
solutions for (N± min , N± max ). Top panel, the 138 datapoints from Paper II; bottom panel, the 69 
datapoints from paper III. 

area that has two peaks with Pks ~ 0.89 connected by a saddle. We don't regard the differences 
within the highest contour to be significant and we adopt the point indicated by a star as the best 
solution; here Pks = 0.86. We adopt its uncertainty as defined by the contours where Pks drops 
to ~ 0.30, which is ~ 32% of the peak value 0.89. This yields 



N ±min = 0.0261K x 10 20 cm -2 (28a) 
N ±max = 2.6t\i x 10 20 cm" 2 (28b) 



These two limits differ by two orders of magnitude! For most of the range of iV_i_ between these 
limits, <p(N±) oc iV^ 1 to a very good approximation. In Figure 7, the smooth solid line in the top 
panel is equation 27 with these values of (N± min , N± max ). The bottom panel shows the cumulative 
distributions, both for the data (in diamonds) and for this derived pdf. To the eye the match looks 
excellent, and this is confirmed by the high value of the K-S probability Pks = 0.86. 
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Figure 7 also shows the distributions for the current magnetically-selected sample 69 data- 
points from Paper III. The top panel shows the histogram in a dotted line. The dotted histogram 
differs significantly from the solid one because it doesn't have the large increase for small Ni os . 
The magnetic selection excludes datapoints whose uncertainty 5B[ os exceeds an upper limit. This 
restriction biases the column densities, because it is impossible to obtain small uncertainties on Bi os 
for small Ni os . Consequently, the pdf of Ni os is cut off at small values. The bottom panel of Figure 
8 shows contour plots of PKs(N± min , N± max ) for this set of 69 datapoints. Clearly, the contours 
peak at a different location. The bias in the magnetically selected points is reflected in the larger 
value of N± min , 0.090 for the magnetically-selected sample versus 0.026 for the unbiased sample. 
The values of N± max are comparable: 2.0 and 2.6 for the small and large samples, respectively. 
The maximum value of Pks = 0.58, indicating good agreement between the data and the analytic 
representation. 

4.3. Derivation of the intrinsic 4>(Vt ur b) from the histogram of observed Vt ur b,los 

The top panel of Figure 9 exhibits the cumulative distribution of the 138 Vt ur b,los from Paper 
II (solid line), together with an adopted cumulative distribution (dashed line). We obtained the fit 
by trial and error: for trial functional forms for (f)(Vt ur b), we used equation 15 to calculate the corre- 
sponding ip(Vturb,ios)- Trial functional forms included (p(Vturb) x n e~ x , with n = (0.25,0.5, l.,2.), 
where x = Y turb . The value n = 1 is significantly better than the others, resulting in Pks = 0.91 
for Vturbfl = 1-06 km s _1 . This functional form is not unique: the function (f>(V tur b) oc x l / 2 e~ x fits 
almost as well, with Pks = 0.81 for V tur b,o = 1-7 km s _1 . These line widths are dispersions. The 
middle panel shows the data histogram and corresponding pdf i/)(Vt ur b,ios), and the arrow indicates 
the median of the observed distribution at V tur b,i os ,i/2 = 1-2 km s -1 . The bottom panel shows 

(f>( V turb)- 

The above paragraph assumes the model in which Vturb is perpendicular to Btot- If we drop 
that assumption, then the line-of-sight distribution of Vt ur b,ios has the same functional form as that 
of Vturb (but is narrower by 3 1 / 2 ). A good fit for the observed distribution is ip(Vt ur b,ios) °c x 2 e~ x , 
i.e. n = 2 in the above paragraph, with Vt ur b,ios,o = 0.45 km s _1 . This line width is dispersion. 

5. DERIVATION OF PREDICTED BIVARIATE OBSERVED DISTRIBUTIONS 

AND COMPARISON WITH DATA 

In this section, we use the previous determinations of the pdf of intrinsic quantities 4>(B to t), 
4>(N±), and 4>(Vt U rb) to determine the expected bivariate distributions of observed pairs ip(B b s , Ni os ), 
V '(Vturb, losi Bobs)-, and tp(V tU rb,iosi ^los)- If the observations were good enough, this would enable us 
to use the correlations among observed quantities to infer astrophysical model parameters. Unfor- 
tunately, we will find that we cannot. 
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Fig. 9. — The top panel shows the fitted and actual cumulative distributions of Vt U rb,los (i n km 
s _1 ). The middle panel shows the corresponding pdfs, i.e. the data histogram and ^(Vt ur b ios)- The 
third panel shows our fitted 4>(V tur b)- The arrow shows the median of the observed histograms. 
Line widths are dispersions. 

5.1. <p(B tot ,N ± ) and *p(B obs , N tos ) 

Here we employ the observationally-derived univariate distributions for Btot and -/Vj_ to pre- 
dict the bivariate distribution of the observed quantities ip(Bi os , Ni os ), for the two cases of B to t 
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perpendicular and parallel to the sheets. However, this single step does not provide a firm basis for 
comparison with the data because the observational uncertainties on Bi os are so large. 

As we saw in §4.1.1, the noise precludes our deriving even the univariate 4>(B tot ), so for discus- 
sion purposes we adopt the EXP FCN of equation 8 with B$ = 3.9 /iG (which provides the required 
median value 6.0 //G) (§4.1.1). Even with this, however, there is additional step: we must predict 
the bivariate distribution of the measured B b s (which includes observational error), i.e. we must 
obtain i/j(B i, s , Ni os ) instead of ip(Bi os , Ni os ). We accomplish this by convolving the conditional 
distribution ip(Bi os \ Ni os ) with a Gaussian having the rms measurement dispersion 3.4 /j,G (§4.1.2), 
and doing this as a function of Ni os . 



5.1.1. Case of B to t Perpendicular to the Sheet 

Here we use equation 17 to predict the observed ip{Bi os ,Ni os ) for the perpendicular case. As 
explained above, for cf)(B tot ) we use the EXP FCN of equation 8 with Bq = 3.9 fiG ; and (f)N ± comes 
from §4.2. The analytic solution is somewhat cumbersome and yields 



tp (Bios, N ios ) 



-[B? M JV? M /(2BgJV ± ^J] _ e -[BlNlJ{2BlN^ mm )] 



Nios < N ±1 
Nos > N ±r 



(29) 



This produces contours in the (Bi os , Ni os ) plane; Figure 10, top panel, shows the results. The 
version including measurement errors is in the middle panel. 

The general trend in the top panel is clear from the discussion of §2: large measured column 
densities Ni os produce small measured fields Bi os because the field lines are perpendicular to the 
sheets. Unfortunately, including the observational uncertainties, as we do in the middle panel, 
obscures this trend. We defer further discussion to §6. 



5.1.2. Case of B tot Parallel to the Sheet 

Here we proceed as in §5.1.1, but use equation 19 instead of equation 17 to predict the ob- 
served ip (Bios, Nios) f° r the parallel case. The math is complicated, so we proceed by calculating 
^(Bios^Nios) using a Monte Carlo simulation. Figure 11 shows the results. The version including 
measurement errors is in the middle panel. 

The general trend in the top panel is clear from the discussion of §2: for Ni os 3> N± max , we 
see sheets more nearly edge-on where Bi os is usually large and the conditional pdf ip(Bi Q s \ Ni os ) 
becomes independent of Ni os - We defer further discussion to §6. 
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Fig. 10. — Top panel, contours from equation 29 of tp(Bi os , Ni os ) for the perpendicular case derived 
using the observationally-derived cf)(B tot , N±). Units are fiG and 10 20 cm~ 2 . The middle panel 
shows ip(B b s , Nios), which includes measurement errors on Bi os . The square marks the univariate 
medians. Bottom panel, the data. 
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Fig. 11. — Top panel, contours of ip(Bi os , Ni os ) from a Monte Carlo simulation for the parallel 
case derived using the observationally-derived (f)(B tot , N±). Units are fiG and 10 20 cm -2 . The 
middle panel shows i(;(B bs, Ni os ), which includes measurement errors on Bi os . The square marks 
the univariate medians. Bottom panel, the 69 datapoints from Paper III. 
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5.2. (t>(V tU rb, B tot ) and ip(V tU rb,ios, Bios) 

As noted in §3.3.1, the bivariate distribution does not depend upon any model regarding the 
geometrical shape of the HI clouds. Figure 12, top panel, shows contours of i){Vt ur b,ios, Bi os ) from a 
Monte Carlo simulation using the observationally-derived <f>(V tur b, B tot ); these Bi os do not include 
measurement errors. The middle panel shows ip(Vt ur b,ios, B b s ) and the bottom panel shows the 
data. Contours are spaced by a factor of two. The eyeball says that the fit to these 69 datapoints is 
not bad, which makes it tempting to conclude that the turbulence is, in fact, perpendicular to the 
magnetic field. Unfortunately, however, the eyeball also says that the distribution of points doesn't 
differ much from one in which \B b s \ and Vt U rb,los are uncorrelated, meaning that such a conclusion 
is unwarranted. 

5.3. 4>{V tur b, N ± ) and ip(V turb ,ios, Nios) 

Here, as above with the relationship between Bi os and Ni os , the predicted bivariate distribution 
depends on whether the magnetic field is perpendicular or parallel to the sheet. Figure 13, top 
panel, shows contours of ip(Vt ur b,ios , Bi os ) for the perpendicular model; the middle panel shows the 
parallel model. The contours are derived from Monte Carlo simulations using the observationally- 
derived 4>(Vt ur b, N±), assuming that two quantities are independent. The bottom panel shows the 
data. 

Unfortunately, the difference between the contours in the top and middle panel is not large 
and it doesn't take a quantitative analysis to tell that the data cannot distinguish between the two 
models. To the eye, the data look in reasonable agreement with both models. 

6. WHICH SHEET MODEL FITS BETTER? 

The key to distinguishing between the two sheet models is the bivariate distributions because 
the univariate distributions don't depend on whether the perpendicular or parallel model reigns. 
Of the three bivariate distributions, two are relevant. 

Figure 13 displays ip(Vt ur b,ios, Ni os ). Unfortunately, the differences between the contours in 
the top panel (perpendicular model) and middle panel (parallel model) are not large. We can 
make a quantitative comparison of these bivariate distributions with the data by using the two- 
dimensional generalization of the K-S test described by Press et al (1997), which provides Pks m a 
similar fashion as does the usual one-dimensional K-S test. The values are Pks = 0.032 and 0.016 
for the perpendicular and parallel cases, respectively. These values are both small, but not so so 
small as to rule out agreement of the data with the model. They differ by a factor of two, but this 
is nowhere near enough of a difference to distinguish between the two cases. Moreover, our basic 
assumption in deriving these distributions is the anisotropy of the turbulent velocity V tur -b, which 



-31 - 



MODEL. IB T 




<>.<P o o 



1 



Fig. 12. — Top panel, contours of ip(Vt ur b,ios, Bi os ) from a Monte Carlo simulation using the 
observationally-derived 4>(Vt ur b, Btot)- Units are km s^ 1 and fiG. The middle panel shows 
ip(Vturb,iosi B bs), which includes measurement errors on Bi os . Bottom panel, the data. 
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Fig. 13. — Contours of ip{Vturb,ios^ios) from a Monte Carlo simulation using the observationally- 
derived 4>(Vt U rb) and (j)(Bt t)- Units are km s _1 and 10 20 cm -2 . Top and middle panels show the 
perpendicular and parallel sheet models, respectively. Bottom panel shows the data. 
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might not be correct. Unfortunately, this bivariate distribution cannot distinguish between the two 
models. 

The bivariate pair tp(Ni os ,Bi os ) is more directly related to the two models because there is 
no additional assumption involving anisotropy of the turbulence. Figures 10 and 11 exhibit the 
pair for the perpendicular and parallel cases, respectively. After including the measurement errors 
(middle panels), the two bivariate distributions look similar, as do the runs of the median B obs l / 2 
versus Ni os . Neither represents the data very well. We can be quantitative and perform the two- 
dimensional K-S test, as above; the perpendicular and parallel cases yield Pks = 0.022 and 0.024, 
respectively. Again, the values are small but not so small as to rule out the models; and, again, 
the test cannot distinguish between the two cases. 

One can argue that the single reliably-detected point at (Ni os ,Bi os ) ~ (7,8 /uG) is more con- 
sistent with the parallel model, for which the two quantities tend to be correlated. This particular 
Gaussian component of 3C142.1 (Paper III) has among the largest values for both quantities, which 
means that for the perpendicular model we need either a very thick sheet and/or a very high in- 
trinsic field strength. But we cannot extend this argument to other sources, whose values are more 
representative. 

A confident choice of model requires more and better data. Unfortunately, this is not likely 
given current instrumentation, because we have already used ~ 800 hours of Arecibo time for this 
project. 

7. ASTROPHYSICAL DISCUSSION: MAGNETIC FIELDS 

7.1. Observational issues 

As discussed in §3.1.1, we observe the line-of-sight component of magnetic field Bi os , which 
is always less than the actual magnetic field B tot . Figure 1 and its associated discussion shows 
that, in practice, we cannot determine the intrinsic distribution (p(B to t) from observations of Bi os . 
Fortunately, the median is well-determined, with \B tottl / 2 \ = 6.0 ± 1.8 (§4.1.1). 

Figure 1 and Table 1 show that the median line-of-sight field strength Bi osl j 2 is less than 
half that of the total field B totl / 2 , with most of the histograms of ip(\Bi os \) peaking at zero. 
This functional behavior of Bi os , which results purely from geometry, is responsible for the large 
number of nondetections in our sample. Moreover, it explains particular cases where Bi os is small. 
A spectacular example is the local-arm field seen against Cas A, Bi os = —0.3 ± 0.6 /xG. This 
surprisingly small result is perfectly consistent with statistical expectation. Of course, we cannot 
rule out that the field actually is really small in any particular case like this, but one needs additional 
data to draw such a conclusion! 

Now consider the large set of magnetic fields observed in 21-cm line emission in morphologically 
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obvious structures, reviewed by Heiles & Crutcher (2005). The term "morphologically obvious" 
means filaments or edge-on sheets. Edge-on sheets should be edge-on shocks in which the field is 
parallel to the sheet. Here, the statistics reverse and favor relatively large Bi os . As explained in 
§3.2.3, as the line of sight becomes parallel to the sheet — i.e., for a morphologically obvious sheet — 
the median B totl / 2 — > 0.7lB tot . For these structures, measured fields are strong, ranging from ~ 5 
to ~ 10 /j,G. This is not inconsistent with a uniform B tot ~ 10 /iG, which is almost a factor of two 
above the median CNM field strength. This suggests that shocks enhance the field strength, but 
not by large factors. 

Finally, compare the CNM median field B to t\/2 with other estimates of field strength. Beck 
(2003) reviews the most recent estimate of field strength derived from synchrotron emission, mini- 
mum energy arguments, measured cosmic ray flux, and polarization. He finds the regular component 
to be ~ 4 fj,G and the total component to be ~ 6 fiG. Pulsars give a much smaller value for the 
regular component, but they provide an underestimate if field and electrons are uncorrelated, as 
is likely (Beck et al 2003); nevertheless, they give about 5 /j,G for the total component (review by 
Heiles 1996). 

The difference between the regular and total components is the fluctuating component, whose 
turbulent spectrum covers a wide range of scales ranging up to at least tens of parsecs (e.g., the 
North Polar Spur). Our CNM structure sizes are typically of order tenths of a parsec (Paper III), 
smaller than much of the magnetic field's turbulence scale range. For this reason we think that it 
is more appropriate to compare the CNM field strengths with the total component given by Beck, 
not the regular one. Our CNM median of ~ 6.0 //G is close to the local Galactic total component 
of ~ 6 fiG. 

7.2. Astrophysical issues 

We find that B tot in the CNM is comparable to the field strength in other ISM components. 
This is at first surprising because the volume density n(HI) in the CNM greatly exceeds that in 
all other interstellar structures except molecular clouds. Flux freezing applies almost rigorously in 
the diffuse gas, even in the HI, and as the interstellar HI changes from CNM to WNM and back 
again, whether by thermal instability or dynamical processes, the transition must occur under the 
constraints imposed by flux freezing. Under the usual flux-freezing ideas, magnetic field strength 
should increase with volume density. If this increase would actually occur, then we would expect 
higher field strengths in the CNM than in other diffuse gas phases because the ISM should exhibit 
approximate thermal pressure equality among the phases. This evidently doesn't happen. In fact, 
this absence of field strength increase for small n(HI) is well known from past studies (e.g. Crutcher, 
Heiles, & Troland 2003, section 3.4), so this is hardly news. 

The field is strong enough to dominate the gas pressure, and therefore the dynamics. With 
a CNM median field B to t ~ 6.0 /j,G, which also applies elsewhere in the interstellar volume, the 
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mag netic pressure is ~ 10400 cm K. This dominates the CNM pressure Pcnm ~ 3000 cm 
(Jenkins & Tripp, 2001; Wolfire et al, 2003). When the field dominates the pressure, it is much less 
affected by the gas pressure or thermodynamic state. 

CNM structures are magnetically subcritical. When gravitation is important, the distinction 
between magnetically sub- and supercritical clouds occurs at ~ 0-38 (Nakano & Nakamura 
1978). In subcritical clouds, the magnetic field dominates gravity and prevents collapse. Our CNM 
sheets have iV|_ < 2.6 and the median B totl / 2 ~ 6.0 fiG, yielding a minimum ratio ~ 2.3, far 
above the supercritical upper limit. Gravity is far from important in these clouds, but if several 
clouds were to coalesce into a gravitationally-important one then magnetic forces would prevent 
gravitational collapse — unless the field were destroyed in the process, e.g. by the annihilation of 
oppositely-directed fields during coalescence of individual clouds. In this sense, the field is strong 
and must dominate the act of star formation in denser clouds that form from less dense interstellar 
gas. 

8. ASTROPHYSICAL DISCUSSION: COLUMN DENSITIES 

It is well accepted that interstellar HI often lies in sheets. Paper II showed this convincingly 
for the CNM structures. In the present paper, we noted that the observed histogram of observed 
column density iVj os falls monotonically with behavior close to iVj os -1 . We then derived the pdf of 
the intrinsic column density for the sheets N± and found that it follows equation 27 between two 
limits, which are rather well defined. Equation 27 behaves much like a iVi" 1 distribution. Figure 
7 convincingly shows that this is a good description of the observations. 

Suppose, first, that the CNM results from shocks. We have in mind the McKee & Ostriker 
(1977) model, in which a supernova shock adiabatically compresses and heats the ambient gas; as 
the gas cools, it does so under roughly constant pressure so its density increases, it cools faster, 
and soon becomes the CNM. As the shocked gas cools it slows. For low-velocity CNM, which is 
primarily what we observe, the swept-up column density depends on the energy injected and the 
ambient density. This dependence is complicated (Cioffi, McKee, & Bertschinger 1988), so we don't 
attempt a detailed discussion. However, the swept-up column densities are not incomparable with 
our upper limit N± max = 2.6 x 10 20 cm~ 2 in equations 27 and 28. Thus, the ~ N±~ 1 dependence 
in equation 27 could be a reflection of the statistical distribution of the relevant function of injected 
energy and ambient density. 

Another possibility is that the CNM arises from kinematical and thermal processes in the 
turbulent interstellar medium. We have in mind structures like those seen in numerical simulations 
of interstellar turbulence such as Vazquez-Semadeni, Gazol, & Scalo (2000; see references in Paper 
II). Consider the "Triad region" discussed in extensively in §8.2 of Paper II and Heiles & Crutcher 
(2005). It has line-of-sight extent ~ 0.05 pc and plane-of-sky extent ~ 20 pc, for an aspect ratio 
~ 200. It also has typical turbulent velocity ~ 1 km s _1 , which makes the line-of-sight crossing 
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time ~ 5 x 10 4 yr. This is very short — interstellar kinematical evolution over human history! If this 
sheet were the result of a slowed shock, it seems remarkable that, in the presence of ISM density 
fluctuations, the distance over which the swept-up column density is accumulated would allow the 
sheet to appear so coherent, and further that it would retain its coherence to be observable. The 
alternative is, we suppose, that CNM structures are transient, and that we can map as extensive 
only those that currently appear to be coherent. 

This N^" 1 behavior does not seem to be a capricious result. Rather, it is a challenge to the 
theorists to reproduce it. Deciding between the above two possibilities, or others, is a matter of 
the explanation reproducing the N^ 1 behavior. 



9. ASTROPHYSICAL DISCUSSION: TURBULENT VELOCITIES 

In §4.3 we modeled the turbulent velocities as being anisotropic, i.e. perpendicular to the 
magnetic field, and derived the intrinsic pdf (ft(Vt ur b) from the observed one ip(Vt ur b,ios)- Here Vturb 
is the the Id component of turbulent velocity; for the 2d case we modeled, the full component is 
2 l l 2 Vturb- In §5.2 we discussed the bivariate distribution ip(Vturb,los, B f, s ). The data fit the model 
with anisotropic turbulence quite well, but the data also fit no correlation quite well, so the results 
are inconclusive. The observations cannot distinguish between anisotropic and isotropic turbulent 
velocities. 

We can, however, discuss topics such as the relative energy densities in turbulent motions and 
magnetism, that is, whether turbulent velocities are superAlfvenic. Doing this requires some care 
in definitions of parameters: we measure line-of-sight turbulent velocity dispersions and these must 
be converted to their 2- or 3-dimensional counterparts; and we must include He as a component of 
the ISM mass density. Finally, we can discuss results in terms of the conventional plasma parameter 
P, in terms of energy densities, or in terms of supersonic and superAlfvenic. 

First we define velocities and the Mach number. Let AVt ur b,id be the Id turbulent velocity 
dispersion and AVth,id be the thermal velocity dispersion. For isotropic turbulence, the full turbu- 
lent velocity is AVt ur b = &Vt ur b,zd = ^Vturb,id- With this isotropic turbulence, the turbulent 
Mach number M tur b is 



MLb = (30) 

Here C s is the velocity of sound; the appropriate sound velocity is the isothermal one because 
thermal equilibrium is reached quickly in the CNM, so 



l2= nkT 
S P 



(31a) 
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The volume density is 

p = (1 + Af He )n H m H (31b) 



where fji e is the fractional abundance of He by number and m# is the mass of the H atom; we 
adopt fue = 0.1. Similarly, for the Alfven velocity Valf, we have 



Vh* = % (32) 



The information propagation velocity perpendicular to the field lines is equal to the Alfven velocity 
Valf- The mean square velocity perpendicular to the field lines is twice the line-of-sight value. 
Consequently, we define the Alfvenic turbulent Mach number MALF,turb as 

2AV 2 

M ALF,turb — 772 \ 66 ) 

ALF 

If MALF,turb > 1, then then shocks will develop; this is the superAlfvenic case. 
Next we define energy densities. The turbulent energy density is 

pFAV t l rb 

murb = 2 \ 6 > 

Here we include the quantity F to allow for anisotropic turbulence. If turbulent velocities are only 
perpendicular to B, then F = 2; isotropic turbulence has F = 3 (the case assumed in paper II). Of 
course, the magnetic energy density is 

Bt t 2 

E ma g = ^ (35) 

The ratio is 

Eturb _ F^Vturb,ld _ F 2 

~~p — 772 ~~ 77 ALF,turb ■ K M ) 

&mag v ALF Z 

Finally, we define the conventional plasma parameter (3th, which compares thermal and mag- 
netic pressures: 

o _ Pth 2 ^ V th,ld 1 + fHe 

Pth = — = -T72 1 + 4f ( 37 ) 

r mag v ALF ' J " e 
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For comparison of turbulent and magnetic effects, we calculate the relevant ratios for the 
following adopted parameter values, which are close to the medians: 



T C nm = 50 K ; (38a) 

= 3000 cm" 3 K ; (38b) 

k 

n{HI) C NM = 54 cm" 3 ; (38c) 

&V turbM = 1.2 km s- 1 ; (38d) 

Btot = 6.0 fiG . (38e) 



Here T is a typical CNM temperature from Paper II. Pcnm is from Jenkins & Tripp (2001) and 
Wolfire et al. (2003). The value for AV tur b,id = 1-2 km s" 1 is the median from §4.3. The value for 
B t ot is the median from §4.1.1. 

These values provide 



M turb = 3.7 ; (39a) 

Pth = 0.29 ; (39b) 

M 2 ALF , turb = 1.3 ; (39c) 

Valf = 1.5 km s- 1 . (39d) 



If the field were small, then turbulence would be isotropic with F = 3. However, it is not so small, 
but neither is it so large that we could say for certain that F = 2. The limits 2 < F < 3 correspond 
to 1.3 < 1^ < 1.9. 

These values should be regarded as representative. Not all CNM clouds have the median 
values, so these parameters have a considerable spread. 

One interesting question is whether any CNM Clouds have high (3th-, i.e. whether any CNM 
clouds have negligible magnetic field. We cannot answer this because @ t h ^ an d neither P nor 
Bf t is attainable for an individual cloud. Moreover, the functional form of Btot is ill determined: 
we cannot even distinguish between all values being the same (a delta-function distribution) and a 
Gaussian (which has a significant population of very high field strengths); see §4.1.1. 
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10. FINAL COMMENTS: TURBULENT AND MAGNETIC ENERGY 

EQUIPARTITION 

Our numbers indicate that magnetism and turbulence are in approximate equipartition. The 
approximate equipartition suggests that turbulence and magnetism are intimately related by mutual 
exchange of energy. In the absence of turbulence, magnetic energies do not dissipate because the 
magnetic field cannot decay on short time scales. But with turbulence, the field may be able to decay 
rapidly (Heitsch & Zweibel 2003; Lazarian & Visniac 1999; Zweibel 2002). Moreover, numerical 
simulations suggest that supersonic turbulence also dissipates rapidly, even when the field is strong 
(MacLow et al 1998). However, it is not obvious to us that these dissipative processes should lead 
to the observed equipartition between turbulence and magnetism. 

We suspect the answer lies in Hennebelle & Perault's (1999) result and Maclow's (private com- 
munication) observation that the CNM components result from the transient nature of turbulent 
flow: the CNM occupies regions where densities are high, produced by converging flows, and the 
density rise is limited by pressure forces. These limiting pressures are magnetic because the gas 
has small [3th, meaning that thermal pressure is negligible and the dynamical equality makes the 
magnetic pressure comparable to the converging ram pressure and leads to apparent equiparition. 

The equipartition looks like a steady-state equilibrium, but it is really a snapshot of time 
varying density fields and our immediate observational view is a statistical result over a large 
sample. In other words, our current observational snapshot shows an ensemble at a given time. 
Against this we compare the numerical simulations, which are stationary in the sense that they 
have been allowed to run long enough that the statistical properties become time-independent. 
Such simulations are also ergodic, with statistical properties over time being equivalent to those 
over space. With this view, the ISM dynamically evolves through turbulence and its properties are 
governed by statistical equilibrium of energy inputs and dissipation. These matters are discussed 
at length in the excellent review by MacLow & Klessen (2004) 

It is a pleasure to acknowledge discussions with Martin Bureau, Jeremy Darling, Avinash Desh- 
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